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Abstract 
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an extension of the Standard Model with 4 generations and the inclusion of four right-handed 
Majorana neutrinos. We obtain a value for the baryon asymmetry of the Universe in accordance 
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1 Introduction 



The attractive mechanism of thermal leptogenesis has been thoroughly investigated in recent 
years. Most of the focus has been on the extended model of the Standard Model (SM) with 3 
right-handed (RH) neutrinos pp. This simple model which can explain the baryon asymmetry 
of the Universe (BAU) and provide appropriate values for neutrino masses and mixing angles 
is based on the out-of-equilibrium decay of a RH neutrino, in a CP violating way, such that 
an asymmetry in the leptonic decay products is produced. The leptonic asymmetry is then 
converted into a baryonic asymmetry due to {B + L)-violating sphaleron interactions which are 
in equilibrium above the electroweak breaking scale. 

From the experimental point of view the recent best fit values of the ratio of the baryon to 
photon density V = tt^ 

r/ = 6.5+0;^ X 10-^0. (1) 

is given by WMAP P]. 

The typical scale at which the standard scenario described above occurs corresponds to 
10^*^ — 10^^ GeV, and the induced neutrino masses arise via the see-saw mechanism (31 1^ |5j 
with a common interaction being at the origin of the asymmetry and mass generation. 

Recently, detailed calculations have been performed solving the corresponding Boltzmann 
equations of the system taking into account the different interactions which contribute both 
to the production and potential washout of the leptonic asymmetry. The implications of these 
analyses have provided strong constraints on the values of the light left-handed neutrino masses 

rrii ~ 0.2 eV [HIIZIIHI and a lower bound on the mass of the lightest RH neutrino M ^ lO^GeV 
ini m ^1 • The more recent work of ref. has included the effect of scattering with gauge 
bosons which had been neglected before and in ref. J2] a comprehensive calculation including 
effects from renormalization group (RG) running, finite temperature background effects, other 
possible decays producing a leptonic asymmetry (decay of a scalar particle), a correction to the 
appropriate subtraction that must be performed for scattering via N-exchange was presented. 
A thorough recalculation of this last issue was done in ref. US]) coinciding with ref. ^2] that 
the washout contribution from this type of diagram had been overestimated. 

The less standard scenario of leptogenesis occurring at the TeV scale has been investigated 
for example in references ^1 El El El iS* JO, 21 . In particular, in ref. ^l] a TeV scale 
model which is a simple extension of the Fukugita-Yanagida (FY) model including one more 
generation to the SM and four gauge singlets was discussed and estimates of the produced 
CP-asymmetry and final baryon asymmetry were obtained. Both supersymmetric and non- 
supersymmetric versions of the model were analyzed. 

In this paper we focus on this model of ref. and solve the corresponding Boltzmann 
equations (BE) in order to reliably obtain the final values of the asymmetry and determine the 
region of parameter space which is available for the model. 

We dedicate section 2 to a careful discussion of the way we set up our calculations of the 
Boltzmann equations making emphasis on a few subtle points. In section 3, we recall the main 
features of the model and establish our notation. Section 4 presents the Boltzmann equations 
valid for this model. Section 5 is devoted to our results and conclusions. 
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2 Preliminary Considerations 



In reference [221 it was pointed out that the relevant object that should be studied in thermal 
leptogenesis is the density matrix p |2S|- In the FY model p is a (3 x 3) matrix with all 
entries different from zero for temperatures greater than lO^^GeV. However, as the temperature 
decreases and each of the different interactions defined by a specific charged Yukawa coupling 
enters equilibrium the corresponding off-diagonal elements of the density matrix go to zero. In 
our model at the TeV scale all of the charged Yukawa couplings are in equilibrium and thus p 
will be a diagonal matrix. The matrix p is normalized such that, 

Trp = Y,YL^^Y.iY,^-^l.) ■ (2) 

a a 

where the sum over a indicates a sum over flavour. 

Another important issue is the relationship between the produced lepton asymmetry and 
the induced B — L asymmetry which is conserved by sphaleron interactions. Again a careful 
consideration of the temperatures and which interactions are in equilibrium changes the equa- 
tion relating these asymmetries, see reference |22j . In general the procedure that is employed is 
through the establishment of equations for the chemical potentials reflecting which interactions 
are in equilibrium. Thus the chemical equilibrium equations allow us to write the asymmetries 
in number densities in terms of some specific asymmetry, say Y^^ . For the Fukugita-Yanagida 
model then generically there would be a (3 x 3) matrix Aa/s of the form j22j 

yB/3-L^=A^f,YL„ (3) 

and the numerical values for the elements of A are determined from the equations for the 
chemical potentials. If all interactions related to charged Yukawa couplings are in equilibrium, 
as well as sphaleron interactions, the equation for the case of N generations simplifies to, 

^-^-" = ("W + 54A^)S*''.-3>'^.. (4) 

or adding over all generations 

Some comments are in order. In the standard scenario of the FY model with three RH 
neutrinos, where the production of the asymmetry occurs at some high scale for which not all 
charged Yukawa couplings are in equilibrium, one should use eq. © to relate the asymmetries. 
The point is that usually what is solved in the system of coupled Boltzmann equations, under 
certain assumptions, is an equation for Yl = J2a ^La-, see below for a discussion on this point, 
and then the asymmetry in — L is calculated using eq. © . From our previous discussion we 
see that it is not quite correct to do this unless some approximations are done concerning the 
different Yl^. 
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For our model at the TeV scale , all charged Yukawa couplings are in equilibrium and it 
would be correct to use eqns. (@J) and ©• Another feature of our model is that the Yukawa 
coupling of the heavy LH neutrino is of the order of 1 to ensure that the masses of the fourth 
generation leptons are heavy enough, and also to give an enhancement to the value of the CP- 
asymmetry produced in the model. However, then the scattering interaction which is associated 
to the Yukawa coupling of the heavy LH neutrino will also be in equilibrium in our model. This 
changes the chemical equilibrium equations such that the relationship between B — L [B) and 
L is modified to be: 

__/l3 8 (13-3jV) \ / 22jV + 13 

Thus, as was estimated in reference Jlj there is no net asymmetry produced in the heavy 
left-handed lepton. Consequently, we can put Y^^ to zero and make the usual assumptions 
about the other Y^^ when constructing the Boltzmann equation for Yl provided we use eq. © 
to relate the produced leptonic asymmetry to the final net B — L asymmetry. Below, and from 
now on, Yl denotes the sum over the asymmetries of only the three light leptonic flavours. 

In reference |12j the effect of renormalization group (RG) running for masses and couplings 
was included. In particular, one of the important effects is that the diagrams with AL = 1 
proportional to gauge couplings can become sizeable compared to those proportional to the 
top Yukawa coupling when considering values for these couplings at the high scale at which 
the lepton asymmetry is being produced. In our model at the TeV scale, we have checked that 
this is a very small effect and we will not include RG effects in our analysis. We also have not 
included the AL = 1 proportional to gauge couplings as we are not in the resonant leptogenesis 
case of ref. [TT] . 

Also in ref. jl2j thermal background effects have been included into the calculation, it was 
pointed out that numerically these effects are particularly relevant when correcting the propa- 
gator of the Higgs scalar field. We include only this thermal correction into our calculation. 

Another important point that should be mentioned, which is crucial for our model, is the 
fact that the actual decay temperature at which the lepton asymmetry is produced is just above 
the electroweak scale; thus one must be careful and check that for our choice of parameters 
the sphaleron interactions are still in equilibrium. Although, in the SM there is no first order 
phase transition, just a crossover, and this is not modified for our model with four generations, 
sphaleron interactions are switched off for temperatures below ~ 100 GeV ,24 . A careful 
calculation of the chemical equilibrium equations in the region close to the electroweak phase 
transition was performed in |25j . 

3 The Model 

The relevant part of the Lagrangian of the model we consider is given by, 

L = LsM + i^R^i^ ^R, - ^{i^R^-^R. + h.c.) - {Xi^L^iPR^<t> + h.c.) , (7) 
^In this model, we denote by a the flavour of the fourth leptonic generation. 



) E (6) 
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where ipRi ^-re two-component spinors describing the right-handed neutrinos and we define a 
Majorana 4-component spinor, Ni = ipR^+ipji. ■ Our index i runs from 1 to 4, and a = e, fi, r, a. 
The a component of corresponds to a left-handed lepton doublet which must satisfy the 
LEP constraints from the Z- width on a fourth left-handed neutrino |2S]- The AJ'^ are Yukawa 
couplings and the field <j) is the SM Higgs boson doublet whose vacuum expectation value is 
denoted by v. 

We work in the basis in which the mass matrix for the right-handed neutrinos M is diagonal 
and real, 

M = diag{Mi,M2,M3, M4) (8) 

and define uid = A*^?;. 

To leading order ^, the induced see-saw neutrino mass matrix for the left-handed neutrinos 
is given by, 

= X'^^M-^X'^v'^ . (9) 

which is diagonalized by the well known PMNS matrix. 

We consider the out-of-equilibrium decay of the lightest of the gauge singlets Nj , which we 
take to be A''i. The decay rate at tree- level is given by 

(AU)^.. '(Kj^aj) rhf^M.M, 

Tn. = i^M,- = V ^ ^M. = . — 10 

Stt ^ ^ Stt ^ 8ttv^ ' ^ ' 



a=l 



where A = A"^ and we have defined 
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{4G) _ 2 



y ^ '^ v\ (11) 



a=l 



r 



To ensure an out-of-equilibrium decay of A^'i it is necessary that ^^^.^j^^^s^ <C 1, where H is the 
Hubble expansion rate at T = Mi. This condition will impose an upper bound on the usual 
effective mass parameter rhi which is defined in for example |27j . 

The CP asymmetry calculated from the interference of the tree diagrams with the one-loop 
diagrams (self-energy and vertex corrections) is |2H] 



/ yim[A'-tA,1^ 
(8^) [AtA.],, V 




(12) 



^Given the fact that in our model, at least one of the RH neutrino masses take values ~ TeV, one has to 
go to next order in the general see-saw formula. However, these corrections are small and not relevant for our 
purpose. 

^To simplify our notation, we write the sums over a = 1 to 4 denoting the sums over all flavours e, fi, r, a, and 

sums over a = 1 to 3 to denote sums over the three lightest flavours. For instance, rn'^'^'' — X]a=i ^"^°a/^°^'^ ^'^ i^ 
the usual effective mass parameter mi defined in the literature for j = 1. 
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where 



fix) = V^[l-(l + x)lnl±^], 



9{x) = (13) 
1 — X 

In reference |14| it was shown that the upper bound of the CP-asymmetry produced in the 
decay of the hghtest right-handed neutrino A'^i is 



3 Miruj 



where denotes the largest eigenvalue of the left-handed neutrino mass matrix m^. Due to 
experimental constraints we require that > 45 GeV, which implies that the bound on ei 
is irrelevant. In general there will be regions of parameter space in which the produced CP 
asymmetry can be very large, thus the final allowed region of parameter space could include 
areas in which the washout processes can be very large as well. In reference ^3] possible 
textures for the y matrix were given to illustrate how to obtain masses consistent with low 
energy data and the appropriate amount of CP asymmetry from the decay of A^i. 



4 Boltzmann Equations 

We will now write the Boltzmann equations for the RH neutrino abundance and the lepton 
densities. With all the above arguments in consideration, the main processes we consider in 
the thermal bath of the early universe are decays, inverse decays of the RH neutrinos, and the 
lepton number violation AL = 1 and AL = 2, Higgs and RH neutrinos exchange scattering 
processes, respectively. 

In our analysis we stick to the case where the asymmetry is due only to the decay of the 
lightest RH neutrino A'^i ^. The first BE, which corresponds to the evolution of the abundance 
of the lightest RH neutrino Y^-^ involving the decays, inverse decays and AL = 1 processes is 
given by 

where z = ^ , s is the entropy density and 7^ . , 7s . are the interaction rates for the decay and 
AL = 1 scattering contributions, respectively. 
The BE for the lepton asymmetry is given by 



dYL z . .Yn, , Yl 



dz sH{Mi) ^ + Yl'' ^ 



(16) 



where ei is the CP violation parameter given by eq. (|12|) and is a function of 7^, and 7^. 
and AL = 2 interaction rate processes, called the washout factor which is responsible for the 

^We have checked that the asymmetry produced by the decay of the second hghtest RH neutrino is very tiny 
in our model compared to the one produced by the decay of A^'i . 
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damping of the produced asymmetry. In eqs. H15() and ()16l) . Y^'^ is the equihbrium number 
density of a particle i, which has a mass rm, given by 



z) 



Ml J ' 



(17) 



where gi is the internal degree of freedom of the particle (5^^, = 2, g^ = 4) and g^, ~ 108. 
Explicit expressions of the interaction rates 7„ , 7„ and are given below. 
The reaction density 7^ is related to the tree level total decay rate of the RH neutrino A'^- of 
eq. (Uni) by 



eg 



7i,, = n^. 



Kiiz) 



(18) 



where Kn{z) are the K type Bessel functions. 

We define the reaction density of any process a + 6—>c + d by 



7' 



dx 



(19) 



where a'^^\x) are the reduced cross sections for the different processes which contribute to the 
Boltzmann equations. For the AL = 1 Higgs boson exchange processes, we have 
4 



0=1 
4 



X — Oo 



X — Oo 



' = 3a„ 



(20) 



where 



3a„ 



a=l 

rrij Ml I X — 



X — 2a,' + 2a/i a,- — 2ah , f x — aj + 
_l — J / 



X — aj + aj^ X — a,-. 



ah 



X 



X — 2a,- + 2a/i a,- — 2a/i ^ / x — aj + ah 
1 — m ' 



X — Qj + Qh X — Uj 



ah 



Tr{XiXu) mj 



/Mi 



- 



\MiJ ' 



ah 



Ml 



(21) 



(22) 



the infrared cutoff ^ is chosen to be equal to 800 GeV due to phenomenological considerations. 
So we have 

ls,=Hf+Hf- (23) 

The expressions of reduced cross sections of AL = 2 RH neutrino exchange processes, which 
washout a part of the asymmetry are given by 



a 



0=1 /3=1 j=l 
4 3 4 



a 



OL=l 13=1 n<j 

4 3 4 



(2) 



a=l (3=1 j=l 



a=l (3=1 n<j 
3 = 1 
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where 



A2) 



1 

1 

2^ 



1 + 



X 



+ 



1 + 



X + a. 



In 



X + a 



a,- 



2tt 



X + aj X + 2a j 
1 1 

+ + 



+ 



On 
1 + — 



2 



J_ \ ( x + an 



1 



In 



a; + 



(2) 
nj 



27r 



1 



X + On + ttj 



In 



(x + aj)(x + an) 

(Xjdn 



+ 



In 



an{x + gj^ 
aj(x + an] 



and 



(x 



o=l /3=1 



then we have 



7x, 



E 



+ 



(2) 

j 



7d, 



(26) 
(27) 

(28) 
, (29) 

(30) 

(31) 



There are several important comments to make at this point. First of ah, notice that 
in eqs. l\W\i . (|2U() and (PT|) . we have summed over all states £a- However in eqs. (OH) and 
H25p terms which correspond to /3 = 4, (that is the a flavour, see figure do not contribute 
to equation ((TC)) because they are multiplied by Yi^ = 0. Indeed, as was mentioned before, 
the processes involving at the same time only this flavour in the initial and final state are in 
thermal equilibrium, and so there is no asymmetry in this flavour. Second, as in the case of FY 
model, one can parametrize 7^ and 7^,. by Mi and rh^^'^'' (see eqs. ((TU)) . (^0)) . (|^T|) 'l. Third, 
the AL = 2 processes can be divided and treated in two separate regimes. The contribution 

in the temperature range from z = 1 to z ~ 10 is proportional to \J m^^^^frS^^'^ instead of 
fhi. For the range of temperature {z >> 1), the dominant contribution at leading order 
would have been proportional to fh? = tr{rnlm^), had the summation on the index /3 gone 
from 1 to 4. In addition, we emphasize that the so-called RIS (real intermediate states) in the 
AL = 2 interactions have to be carefully subtracted, to avoid double counting in the Boltzmann 
equations. In our calculation, we have followed ref. JH] where the authors have shown how 



the appropriate (- 



'87l3, 



in eq. (|31j) ) subtraction must be done. 



5 Results and Conclusions 

5.1 Constraints on the fourth generation 

Our model is based on the addition of a fourth doublet of leptons and quarks (in order to 
be anomaly-free) to the particle content of the Standard Model, and of four RH neutrino 
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4 iaiD </> 




Figure 1: s and t channel AL = 2 Majorana exchange diagrams. 



singlets. The masses and couplings of the singlets are constrained by the left handed neutrino 
masses (from solar, atmospheric and WMAP data) and the right amount of CP asymmetry 
ei. Concerning the fourth generation of leptons and quarks, the direct constraint is that they 
must be heavier than There are stronger constraints from electroweak precision tests: 

for example, the lower bound on the charged lepton a from LEP II is approximately 80 GeV 
j29j . As is well known, the SM central value of the Higgs mass is lower than the direct lower 
bound set by LEP II [SHj and the existence of a fourth generation is one possible way to 
increase the Higgs mass [311 1321 IH^ . The Z-lineshape versus fourth generation masses have 
been extensively studied in for example ref. [H^ where we can see that the fourth left-handed 
neutrino is excluded at 95% CL if its mass is lower than 46.7 it 0.2 GeV. In our analysis the left 
handed neutrino masses are always constrained in order to fit the data (solar + atmospheric 
+ WMAP) as well as this latter bound for the 4th LH neutrino. There are other constraints 
on the fourth generation, the generalized CKM matrix elements, etc. which can be found in 
j2f)j and in ref. jSS]) where present and future experimental searches at Tevatron and LHC are 
widely discussed. 



5.2 Numerical solutions of BE 



In our numerical analysis we use one of the possible textures for the Yukawa matrix(Ajy) in 
the four generation case which can produce the needed amount of CP-asymmetry and fit the 
neutrino data, 



K = C 



( 



\ 



\ 



a 



1/C / 



(32) 



where C, e and a are parameters (e and a are complex numbers). This texture will induce to 
first order the following mass matrix for the three lightest LH neutrinos 



VTLy (X 




(33) 



This is a typical texture for a hierarchical spectrum that has to fit the oscillation data (atmo- 
spheric, solar, CHOOZ) and the constraints on the absolute mass. 
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To show the feasibihty of our scenario we chose different values of the parameters C, e, a 
and values of the RH neutrino masses, Mj, i = 1, 4 and we constrain the inputs given the data 
for the LH neutrino masses and the out-of- equilibrium condition. It is precisely a hierarchy 
in the Yukawa couplings which allows for an out-of-equlibrium decay without a suppression of 
the induced CP-asymmetry. 

In figurelSl we illustrate for a given set of input parameters the different thermally averaged 
reaction rates contributing to BE as a function oi z = 

r^ = ^, X = A 5, AL = 2. (34) 

It is clear from this plot that for this set of parameters, and this is true for a wide range of 
parameter space, all rates at z = 1 fulfill the out-of-equilibrium condition (i.e. Vx < H{z = 1)), 
and so the expected washout effect due to the AL = 2 processes will be small. The parameters 
chosen for this illustration are: Mi = 450 GeV, M2 = 2 x 10^ GeV, M3 = 10^ GeV, M4 = 605 
GeV, ei ~ 4.2 x IQ-^ mf*^^ = 2.7 x IQ-^ eV and mf'^^ = 3.6 x 10"^ eV. 

Although a complete scan of the allowed parameter space varying both the values of the 
RH neutrino masses and the Yukawa couplings has not been done we mention that it is not 
completely trivial to satisfy all constraints when we vary the parameters. For example, for 
fixed values of the Yukawa couplings, if we let the second lightest RH neutrino be heavier 
then the heavier LH neutrino becomes lighter than Mz/2. On the other hand, again for fixed 
values of the Yukawa couplings, when the masses of the two heaviest RH neutrinos increase 
or decrease, it becomes difficult to have acceptable light LH neutrino masses. Nevertheless, a 
more detailed anaylsis can provide other spectra for the RH neutrinos which are also acceptable. 

Figure El represents the solution of the BE, abundance and B — L asymmetry, as a function 
of z. The generated value of the baryon asymmetry is 77^ ~ 6 x 10^^*^. Applying the see-saw 
mechanism to our model for the chosen values of the parameters, we obtain a heavy left-handed 
neutrino with a mass above 48 GeV, which is consistent with the fourth generation constraints 
(see section l5?T|) and the three light neutrino masses are of the order of 10~^ eV to a few 10~^ 
eV. 



We have presented the solutions to the coupled system of Boltzmann equations for our TeV 
scale model of thermal leptogenesis. We have carefully considered the effect of the interactions 
involving the heavy fourth generation leptonic fields and consistently written the BE which 
contribute to the final baryon asymmetry together with the appropriate conversion factor from 
the lepton asymmetry to the baryonic one. Our results show that in this simple extension of 
the Standard Model it is possible to produce the right amount of the BAU in a generic way. 
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